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ABSll^CT 


The  turbuL'iice  generated  by  waves  breaking  on  a  natural  beach  is 
examined  using  hotfilm  anemometer  data.  Turbulence  intensity  is  estimated 
from  the  dissipation  rate  and  an  appropriate  length  scale  (a  fraction  of  the 
water  depth).  The  dissipation  rates  are  determined  from  wavenumber  spectra 
found  by  applying  Taylor's  hypothesis  to  frequency  spectra  of  short  (1/8  s) 
hotfilm  time  series.  The  resulting  Froude-scaled  turbulence  intensities  are 
relatively  uniform  throughout  the  water  column  and  are  similar  in  vertical 
structure  but  lower  in  magnitude  than  in  existing  laboratory  studies.  The 
magnitudes  of  the  turbulence  intensity  observed  in  both  the  field  and 
laboratory  are  consistent  with  an  existing  macroscopic  model  of  bore 
dissipation  in  the  surf  zone.  Scaling  by  this  bore  model  relates  turbulence 
intensity  levels  of  monochromatic  waves  in  small-scale  laboratory  experiments 
to  random  waves  in  the  natural  surf  zone. 
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1.  INTRODUCTION 


Shoreward  propagating  sea  and  swell  tnergy  is  transformed  in  the  surf 
zone  into  motions  of  many  different  types  and  scales,  including  steady 
currents,  low-frequency  waves,  organized  vortical  flows  and 
high-ReynoIds-number  turbulence  [Battjes,  1988].  Much  attention  has  been 
devoted  to  the  sea-swell  and  infragravity  frequency  bands  (nominally  0.005  Hz 
to  0.5  Hz)  because  these  motions  dominate  the  sea-surface  elevation  and 
velocity  spectra  in  typical  namral  surf  zones.  Less  is  knowm  about  higher 
frequency  turbulent  motions,  which  contain  only  a  small  fraction  of  the  total 
kinetic  energy  of  the  nearshore  velocity  field,  but  usually  dissipate  most  of  the 
shoreward  energy  flux.  Although  many  laboratory  studies  of  surf  zone 
turbulence  have  been  conducted  using  breaking  progressive  waves  [e.g., 
Nadaoka  &  Kondoh,  1982,  Hattori  and  Aono,  1985,  and  Stive,  1980],  hydraulic 
jumps  [Resch  and  Leutheusser,  1972],  and  solitons  [Skjelbreia,  1987], 
experiments  in  the  natural  surf  zone  have  been  hindered  by  instrumentation 
difficulties. 

New  measurements  of  turbulence  intensity  in  the  natural  surf  zone  are 
described  in  Section  2.  In  Section  3,  the  wavenumber  spectrum,  is 
obtained  from  a  point  measurement  in  this  strongly  oscillatory  flow  using 
Taylor's  hypothesis  (subject  to  a  condition  similar  to  Lin's  [1953]  criterion  for 
the  application  of  Taylor's  hypothesis  in  shear  flow).  From  the  inertial 
subrange  of  <1>(k)  we  estimate  the  dissipation  rate  e  and  then  the  turbulence 
intensity  using  a  form  of  the  classical  relationship  e  =  u'^/l  where  u'  is  the  rms 
turbulence  intensity  and  I  is  the  energy-containing  eddy  scale.  In  Section  4  we 
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show  that  Froude-sca!ed  sail  zone  turbulence  levels  are  smaller  in  the  natural 


surf  zone  than  in  existing  laboratory  studies.  An  existing  model  tor  the 
dissipation  in  a  bore  shows  that  the  differences  occur  because  the  wave 
frequencies,  wave-height  to  water-depth  ratios,  breaking  intensities,  and 
percentages  of  broken  waves  are  substantially  different.  The  turbulence  levels 
observed  in  both  the  laboratory  and  field  are  consistent  with  the  bore  model. 

2.  EXPERIMENTS 

Experiments  were  conducted  at  Scripps  Beach,  La  Jolla,  California 
during  six  days  in  March  and  April  1992  in  water  depths  between  2S  and  274 
cm.  This  is  a  fine-grained  gently  sloping  (about  1  in  40)  sandy  beach  \vith 
relatively  small  alongshore  depth  changes.  Waves  broke  by  both  spilling  and 
plunging.  Significant  waveheights  in  7  m  water  depth  (a  few  hundred  meters 
offshore)  ranged  from  50  to  120  cm. 

Data  were  collected  from  three  vertically  separated  hotfilm 
anemometers,  two  electromagnetic  (EM)  current  meters  with  4  cm  spherical 
probes,  one  pressure  sensor,  and  a  videocamera  recorder  (VCR).  The  hotfilms 
were  Thermo-Systems,  Inc.  (TSI)  model  1755  constant-temperature 
anemometers  with  1210-60VV  cylindrical  quartz-coated  platinum  hotfilm 
probes.  Although  the  most  rugged  of  the  cylindrical  TSI  probes,  they  are 
fragile.  The  more  rugged  conical  1230W  probes  have  less  desirable  symmetry 
and  gain  characteristics.  Relatively  low  operating  temperatures  (40  °C, 
producing  an  overheat  ratio  of  roughly  8%)  were  used  to  avoid  bubble 
formation  on  the  probe. 
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All  instruments  except  the  VCR  were  mounted  on  a  steel  pipe  frame 
(Figure  1)  which  was  lowered  daily  by  crane  from  Scripps  pier  and  anchored 
to  the  sea  bed  10  m  up-drift  (in  the  alongshore  direction)  from  the  pier.  I  he 
hotfilms  were  located  40  cm  up-drift  from  the  EM  current  meters  and  the 
pressure  sensor  was  buried  about  10  cm  in  the  sea  bed.  The  VCR  was  used  tt) 
determine  the  frequency  and  intensity  of  wave  breaking.  The  time 
synchronization  between  the  videocamera  and  the  hotfilms  was  inadequate  to 
determine  which  individual  waves  were  broken.  Anti-aliasing  analog  filters 
were  applied  to  the  hotfilm  and  current  meter  signals  before  digitization. 

There  were  70  data  runs,  each  of  512-second  duration.  From  the 
resulting  potential  210  hotfilm  time  series,  those  with  broken  probes, 
instruments  out  of  the  water,  or  questionable  calibrations  were  rejected, 
leaving  125  usable  hotfilm  time  series.  The  vertical  positions  of  the  hotfilms 
are  shown  in  Figure  2.  While  the  broad  wave  height  distribution  produces  no 
consistent  "break  point",  breaking  waves  were  infrequent  deeper  than  h/H„  = 
3.0. 

Estimates  of  the  aeration  in  the  air  bubble  region  (foam)  on  the  leading 
face  of  a  bore  range  from  a  few  percent  to  as  much  as  40  percent  [Fuhrboter, 
1970].  The  signal  from  a  constant-temperature  hotfilm  anemometer  is 
corrupted  by  the  bubbles  because  air  has  much  lower  heat  capacity  than  water 
The  bubble  contamination  was  removed  from  the  hotfilm  signal,  before 
calibration,  as  described  in  Appendix  A. 

Because  of  the  difficulty  of  maintaining  stable  calibrations  between  the 
calibration  facility  and  the  ocean,  the  hotfilms  were  calibrated  "in  situ"  using 
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EM  current  meters  located  at  the  same  vertical  elevation  as  the  upper  and 
lower  hotfilm  probes  (See  Appendix  B).  The  mean  of  the  EM  current  meter 
signals  was  used  to  calibrate  the  nearly  equidistant  middle  hotfilm.  Figure  3 
shows  time  series  of  sea  surface  elevation  and  calibrated  hotfilm  speeds  for 
cases  in  which  most  (panels  a-d)  and  few  (panels  e-g)  waves  were  broken. 

The  increase  in  high-frequency  velocity  fluctuations  from  outside  to  inside  the 
surf  zone  (compare  panels  g  and  c)  is  much  larger  than  the  corresponding 
increase  in  wave  height  and  orbital  speeds  (compare  e  and  f  with  a  and  b), 
suggesting  that  breaking  waves  (not  vortex  shedding  or  vibration  of  the 
instruments)  generate  the  high  frequency  fluctuations. 

3.  RESULTS 

Methods  of  Estimating  Turbulence  Intensity 

Turbulence  in  steady  free-stream  flow  or  to  wed-body  experiments  is 
often  defined  as  the  fluctuation  about  the  mean  velocity  [e.g.,  Hinze,  1976].  In 
a  laboratory  surf  zone  with  monochromatic  plane  waves,  turbulence  has  been 
defined  as  the  deviation  from  an  ensemble  average  of  velocities  at  the  same 
wave  phase  [Flick,  Guza,  and  Inman,  1981].  The  stochastic  nature  of  natural 
orbital  wave  velocities  precludes  use  of  this  definition. 

Thornton  [1979]  separated  surf  zone  orbital  wave  energy  from 
turbulence  by  defining  the  wave  orbital  motions  as  the  velocity  fluctuations 
coherent  with  the  sea  surface  elevation  and  assuming  that  all  incoherent 
velocities  are  turbulence.  However,  this  definition  does  not  include  as 
turbulence  the  largest  scale  eddies  which  do,  in  fact,  influence  the  sea  surface 
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[Nadaoka,  Hino,  and  Koyano,  1989].  These  large  eddies  may  contribute 
significantly  to  the  Reynold's  stresses.  Further  disadvantages  are  that 
nonlinearity  [e.g.  Flick,  Guza,  and  Inman,  1981]  and  directional  spreading 
[Herbers,  Lowe,  and  Guza,  1991]  in  a  non*turbulent  wave  field  also  reduce  the 
coherence  between  sea  surface  elevation  and  velocity.  While  the  relative 
underestimation  of  the  more  energetic  orbital  flow  is  small,  overestimation  of 
the  relatively  low  turbulent  energy  levels  may  be  large. 

The  hotfilm  frequency  spectra  (Figure  4)  typically  exhibit  distinct 
turbulent  (expected  inertial  subrange  slope  is  f and  orbital  wave  (expected 
spectral  slope  is  f  ^  [Thornton,  1979])  regimes.  However,  the  transition  (where 
the  slope  changes)  between  the  two  regimes  indicates  only  where  dominance 
changes,  not  the  low-frequency  end  of  the  turbulent  inertial  range.  There  is  an 
overlap  of  unknown  extent,  which  prevents  accurate  estimation  of  turbulence 
intensity  by  high-pass  filtering  above  a  single  "cutoff  frequency”  [Nadaoka  and 
Kondoh,  1982].  Because  the  coherence  and  cutoff  methods  are  inaccurate,  we 
calculate  the  turbulence  intensity  from  the  dissipation  rate  as  described  below. 

Taylor's  Hypothesis 

The  dissipation  rate,  e,  of  a  one-dimensional  wavenumber  spectrum 
<I>(k)  may  be  found  from  the  universal  form  of  the  inertial  subrange 

4)(k)= 


(1) 


where  k  is  the  wavenumber  magnitude  and  a  (nominally  0.5)  is  the  one- 
dimensional  Kolmogorov  constant.  We  first  estimate  from  the  measured 
frequency  spectrum  <t>(f). 

Laboratory  experiments  commonly  concern  relatively  weak  turbulence 
intensity  (u')  in  a  spatially  and  temporally  steady  advective  flow  (U); 

U  >  u'.  (2) 


In  this  case  Taylor's  hypothesis  can  be  used  to  convert  the  entire  frequency 
spectrum  (from  time  series  measured  at  a  fixed  point)  to  the  desired 
wavenumber  spectrum: 


<t)(K) 


4)0 

271/1/ 


(3) 


where 


K  - 


27:/ 

~Lr’ 


(4) 


When  applied  in  situations  with  an  unsteady  advective  flow  field,  U  signifies 
all  non-turbulent  flow:  orbital  wave,  lower  frequency,  and  mean.  Because  the 
orbital  wave  velocities  that  dominate  the  nearshore  spectrum  change  over  a 
wave  period,  we  use  1/8  s  sections  of  the  512  s  records  and  choose  only  those 
where  (2)  is  satisfied.  To  verify  (2),  u'  was  crudely  estimated  by  high-pass 
filtering  the  full  512  s  time  series  at  approximately  the  transition  between 
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the  orbital-wave  and  turbulent  inertial  regimes.  Using  the  deimition  t.  =  L  ^ 
{g/h)‘ '  \  (I  ~  h  is  the  energetic  eddy  scale  and  L'  -  (gh)‘  '■  is  a  common  sur! 
zone  scaling  velocity),  U  fell  near  the  slope  break  in  the  velocity  spectrum  (e.g. 
Figure  4). 

Changes  in  the  wave-dominated  flow  field  during  the  1 /H  s  records  wii! 
also  distort  the  f-K  transformation  (4).  Lin  (1953]  considered  a  similar  tvpe  of 
distortion  in  strong  shear  flows  and  suggested  that  Taylor's  hypothesis  may  be 
applied  to  the  high  wavenumber  turbulent  flow  fluctuations  only  if  the 
nonuniformity  of  the  convection  velocity  over  the  eddy  size  is  small  compared 
to  the  convection  velocity  itself: 

iTldU/dz 

K  >  - - — ,  (a) 

U 

where  z  is  the  vertical  (transverse)  axis.  The  temporal  unsteadiness  in  the 
present  case  is  analogous  to  the  spatial  shear.  We  modify  Lin's  criterion  (5), 
for  use  in  unsteady  flow^  to  simply 

f/  >  At/,  (6) 

w^here  AU  is  the  magnitude  of  the  change,  over  the  1  /8  s  record,  of  the  low- 
pass  filtered  cross-shore  velocity  U. 

One  further  criterion,  unrelated  to  Taylor's  hypothesis,  insures  that  only 
the  u  component  of  the  turbulent  velocity  is  used  in  the  dissipation  rate 
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estimate,  (u,  v,  and  w  denote  tlow  components  in  the  cross-shore,  iongslrore, 
and  vertical  directions,  respectively.)  This  criterion  (see  Appendix  C), 

U  >  V,W  (7 

is  necessary  because  the  hotfilm  probe  senses  flow  from  ail  directions.  We 
have  no  independent  measure  of  W,  the  low-frequency  vertical  flow,  but 
expect  that  U  »  W  for  long  waves  in  the  shallow  waters  of  the  nearshore. 

U,  V  and  AU  were  estimated  for  each  1/8  s  data  segment,  and  u'(t)  was 
estimated  for  each  512  s  data  segment  as  described  above.  When  a  1/8  s  data 
segment  did  not  satisfy  all  criteria  (Equations  2,6,  and  7),  the  1/8  s  data 
"window"  was  moved  in  increments  of  1/32  s  through  the  data  until  the  next 
usable  section  was  found.  The  >>  factor  for  all  criteria  was  chosen  as  5. 
Varying  the  factor  from  4  to  7  typically  produced  10%  variation  in  the 
resulting  u'  estimate.  About  45%  of  the  data  passed  all  criteria,  with  extremes 
of  5%i  and  65%  for  individual  512  second  records. 

Dissipation  Rate 

Each  selected  1/8  s  data  segment  was  detrended,  a  frequency  spectrum 
calculated,  and  a  wavenumber  spectrum  found  from  (3)  and  (4).  A  dissipation 
rate  was  calculated  by  applying  (1)  to  the  best  fit  line  through  the 
wavenumber  spectrum.  (Perhaps  because  of  the  short  time  series,  there  were 
wide  range  of  spectral  slopes.  The  mean  for  all  1/8  s  pieces  was  d.25,  less 
than  the  expected  -5/3,  possibly  owing  to  the  effect  of  bubbles  on  the 
dissipation  rate  [Wang,  1985].) 
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A  plausible  next  step  would  be  to  average  the  1/8  second  e  values 
within  each  512  s  record.  However,  e  in  the  natural  surl  zone  is  higlily 
intermittent  (e.g..  Figure  3(c))  and  Oboukhov  [1962]  suggested  that 
intermittency  in  otherwise  homogeneous  turbulence  produces  an 
approximately  lognormal  distribution  of  dissipation  rate.  The  expected  value, 
<e>,  of  a  lognormally  distributed  dissipation  rate  does  not  etqual  the  mean, 
mode,  or  median  £,  as  it  would  for  a  normally  distributed  variable.  If  e  is  a 
lognormally  distributed  random  variable,  e.g.,  Z=In£  is  normally  distributed, 
then 

<e>  =  (g  +  o^ine/2)  (S) 

where  p  and  cTi^e,  are  the  mean  and  variance  of  Z.  Surf  zone  dissipation  rates 
are  generally  lognormally  distributed  and  the  sample  sizes  small,  so  a 
graphical  procedure  [Baker  and  Gibson,  1987,  see  Figure  5]  was  used  to 
estimate  p  and  cr,„5.  and  thus  <e>  for  each  512  s  data  segment.  Uncertainty  in 
the  resulting  <e>  introduced  by  the  estimation  process  produces  uncertainty  in 
the  turbulence  intensities  ranging  from  ±3%  to  ±15%,  with  a  mean  of  ±5%. 
Record  lengths  of  1/16  and  1/4  s  produced  turbulence  intensities  within  ±10b' 
of  the  1/8  s  values. 

Kolmogorov  [1962]  defined  as  the  intermittency.  Natural  surf  zone 
intermittencies  ranged  from  about  2  to  12,  with  most  values  falling  between  3 
and  8,  while  deep  ocean  values  ranged  from  3  to  7  [Baker  and  Gibson,  1987]. 
The  highest  intermittencies  occurred  where  very  few  waves  were  broken. 
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Dissipation  rates  (from  the  graphical  procedure)  for  512  s  records  in  the* 
natural  surf  zone  ranged  from  0.5  to  500  em'/s^.  Dissipation  rates  in  100  cm 
water  depth  were  typically  of  order  100  cm-/s-\  compared  to  10  ^  to  10  "  cm-/s* 
in  the  equatorial  undercurrent  [Crawford  and  Osborne,  1980]  and  10  '  to  1 
cm^/s^  in  the  tidal  channel  observed  by  Grant,  Stewart,  and  Moilliet,  [1962]. 

Figure  6  shows  normalized  dissipation  rates  in  and  outside  of  the  surf 
zone.  Dissipation  rates  within  the  surf  zone  were  generally  much  larger  than 
seaward  of  the  surf  zone.  There  is  a  weak  dependence  of  the  dissipation  rate 
on  h/Hg  within  the  surf  zone. 

Dissipation  rates  also  were  calculated  using  (1)  applied  to  wavenumber 
spectra  produced  by  applying  Taylor's  hypothesis  to  the  frequency  spectra  of 
512  s  runs  with  U  equal  to  the  r.m.s.  orbital  velocity  [L-umley  and  Terray, 
1983].  This  is  a  crude  approximation  because  the  criteria  for  Taylor's 
hypothesis  are  sometimes  grossly  violated.  Dissipation  rates  were  also 
calculated  by  integrating  the  dissipation  spectrum,  formed  from  the  512  s 
wavenumber  spectrum,  over  the  approximate  inertial  and  viscous  subranges 
[Hinze,  1976].  Both  methods  provided  dissipation  rates  within  roughly  a 
factor  of  two  of  the  <e>  values  calculated  from  the  1/8  s  records  with  the 
graphical  procedure. 

Turbulence  Intensity 

Turbulence  intensity,  u',  was  estimated  from  the  dissipation  rate,  <£>, 

using 
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u'  -  (r4)(K)5K)‘'2 
■^*0 


(9) 


where  <1>(k)  is  defined  by  (1),  giving 

u/^(lay/2(A)W  (10) 

2  K(, 


Svendsen  [1987]  suggested  that  the  turbulent  length  scale  1=2k/k,  is  between 
0.2h  and  0.3h.  We  used  /  =  0.25h  (note  the  weak  sensitivity  of  u'  to  the  exact 
value  of  1),  except  for  sensors  closer  to  the  bed  than  0.25h  that  were  influenced 
by  the  bottom  boundary  layer  (i.e.,  <e>  at  the  bottom  sensor  was  larger  than 
<e>  at  the  middle  sensor),  in  which  case  /  was  set  equal  to  the  distance  to  the 
bed. 

For  three-dimensional  turbulence,  the  inertial  dissipation  estimator  used 
is  consistent  with  Gibson's  [1990]  definition  of  turbulence;  the  inertial  vortical 
forces  are  larger  than  the  damping  forces  because  buoyancy  forces  in  the  surf 
zone  are  small  relative  to  the  inertial  vortex  forces.  However,  the  inertial 
dissipation  method  does  not  include  as  turbulence  the  large  two-dimensional 
eddies  and  "eddy-like  flow"  [Nadaoka,  Hino,  and  Koyano,  1989]  because  they 
are  not  dissipative.  These  flows  may  contribute  substantially  to  the  Reynold's 
stresses. 

Figure  7(a)  shows  that  the  normalized  turbulence  intensities  within  the 
surf  zone  (h/Hg  <  3.0)  are  roughly  10%  of  the  orbital  wave  velocities. 
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Turbulence  intensities  are  smallest  near  the  bed  and  slowly  increase  over  the 
water  column,  consistent  with  strong  mixing  oi  turbulence  irom  the  surtace 
downward. 

Turbulence  intensities  were  also  estimated  from  the  coherence  between 
sea-surface  elevation  and  cross-shore  orbital  velocity.  They  were  generally 
greater  than  30%  of  orbital  wave  velocities  (compared  to  10'.7  by  the 
dissipation  method)  and  independent  of  fraction  of  broken  waves.  The 
directional  spread  of  the  sea-swell  wave  field  reduces  the  coherence  and 
severely  biases  the  estimate  of  u'. 

Turbulence  intensities  were  compared  to  a  model,  based  on  the  energv 
balance  in  a  hydraulic  jump,  which  yields  a  depth-averaged  dissipation  rate 
[Thornton  and  Guza,  1983], 

<e>=—pgy^hfB^  (11) 

4 


where  f  is  the  wave  frequency,  y  is  the  wave-height  to  water-depth  ratio,  and  B 
is  a  breaker  coefficient  related  to  the  intensity  of  wave  breaking.  Combining 
(10)  and  (11),  using  k„  =  2jr/(0.25h),  and  following  Thornton  and  Guza  [1983], 

(12) 

o  4k 


where  w^,  is  the  fraction  of  broken  waves.  B  and  w,,  were  determined  from  the 
videocamera  records  for  each  512  s  data  run.  Past  estimates  of  B,  found  by 
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fitting  the  observed  wave  height  decay  across  both  lab  and  field  surf  zones  to 
model  predictions  based  on  (11),  are  about  1.0  [Thornton  and  Ciuza,  10S3j, 
However,  the  precise  relationship  of  B  to  breaking  wa\'e  properties  is 
unknown.  We  estimated  B  as  the  fraction  of  the  bore  face  that  is  foam- 
covered;  the  range  was  roughly  0.4  to  0.9  with  a  mean  value  for  all  data  of  0.7 
The  broken-wave  fractions  varied  from  less  than  0. 1  to  near  1.0  and  the  wave- 
height  to  water-depth  ratio,  y,  varied  from  0.2  to  0.6. 

The  ratio  u'/u\,  is  surprisingly  close  to  1  (Figure  7(b)),  considering  that 
(11)  provides  only  a  rough  estimate  of  the  dissipation  rate.  The  mean 
deviation  (over  the  vertical)  of  \x' /\x\  from  1.0  is  about  -f-.05,  less  than  the 
scatter  in  the  estimates  themselves  at  a  given  vertical  position.  Turbulence 
intensities  were  also  Froude  scaled  (i.e.,  normalized  by  (gh)‘  ’,  [Svendsen, 
1987])  in  Figure  7(c)  for  comparison  with  laboratory  data  in  the  following 
section. 

4.  COMPARISON  TO  LABORATORY  RESULTS 

The  present  results  are  are  compared  to  previous  laboratory  studies  bv 
Stive  [1980]  and  Hattori  and  Aono  [1985]  (H  &  A)  in  Figure  8,  where  u'  is 
Froude  scaled.  While  the  weak  dependence  on  vertical  position  is  similar  in 
the  field  and  laboratory  profiles,  the  scaled  field  turbulence  intensity  is  less 
than  1/2  of  the  laboratory  values. 

These  laboratory  studies  used  monochromatic  and  unidirectional  w<i\-es, 
each  of  which  breaks  with  about  the  same  intensity.  Howe\’er,  nalurallv 
occuring  waves  are  stochastic  with  a  mix  of  unbroken  and  broken  waves  of 
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various  heights  and  frequencies.  These  complexities  are  heuristically  jncludeLi 
in  the  bore-model  based  estimates  of  turbulence  intensity  (Figure  7b).  The 
bore  model  can  also  be  applied  to  monochromatic  waves  by  setting  both  the 
breaker  coefficient  and  broken-wave  fraction  equal  to  1.0  and  using  observed 
values  of  y.  This  scaling  collapses  the  u'  values  from  lab  and  field  to  within  a 
factor  of  two  (Figure  9),  an  improvement  over  Froude  scaling  (Figure  8).  The 
mean  deviation  (over  the  vertical)  of  u'/u'p  from  unity  is  about  +0.35  for 
Stive's  results  and  -0.20  for  Hattori  and  Aono,  compared  to  +0.05  for  the  field 
data,  (Figure  9(a)). 

Stive  separated  turbulent  and  wave  motions  by  ensemble  phase 
averaging  which  probably  slightly  overestimates  the  turbulence  intensity 
because  of  irregularities  unrelated  to  turbulence  in  paddle-generated  waves 
[Svendsen,  1987].  This  is  consistent  with  u'/u'^  >  1  as  observed  for  the  Stive 
data  (Figure  9a). 

The  bore-model  estimates  of  u'^  for  Hattori  and  Aono  may  be  biased 
low  because  their  data  were  taken  over  a  flat  bed  shoreward  of  a  1  /20  beach 
slope  which  extended  from  deep  water  to  near  the  break  point.  In  their 
experiment,  y  decreased  from  roughly  0.8  near  the  break  point  to  0.4  as  the 
bores  progressed  shoreward  (this  variation  was  included  in  (12)).  It  is  possible 
that  B  also  decreased,  but  B  was  assumed  1.0  in  (12).  H  &  A  estimated  u'  with 
sea-surface  elevation  coherence.  In  a  laboratory  wave  channel  where 
directional  spreading  is  probably  unimportant,  the  method  may  underestimate 
u'  by  not  including  the  largest  eddies  as  turbulence.  Overestimates  of  B  and 
underestimates  of  u'  would  both  tend  to  give  u'/u'^  <  1  as  observ^ed  (Figure 
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9b).  Nadaokn  and  Kondoh's  [1982]  results  (nut  shown)  were  biased  vers'  low. 
(about  30%  of  VI 'p),  consistent  with  their  expectations  that  the  high-pass  filter 
method  used  to  separate  turbulence  from  orbital  velocities  introduces  large 
errors. 

5.  CONCLUSIONS 

Turbulence  intensities  in  the  natural  surf  zone  were  estimated  using 
dissipation  rates  and  a  characteristic  eddy  length  scale  equal  to  1/4  of  the 
water  depth.  The  dissipation  rates  were  determined  from  wavenumber  spectra 
found  by  applying  Taylor's  hypothesis  to  frequency  spectra  of  1/8  s  hotfilm 
time  series.  The  dissipation  rates  were  intermittent,  requiring  lognormal 
statistics  to  determine  expected  values.  The  measuiements  show  that  surf 
zone  turbulence  intensities  are  smallest  near  the  sea  floor  and  slowly  increase 
over  the  water  column,  consistent  with  strong  mixing  from  the  surface 
downward.  The  turbulence  intensities  agree  very  well  with  predictions  from  a 
macroscopic  bore  dissipation  model  which  includes  the  effects  v»f  variations  in 
wave  frequency,  broken-wave  fraction,  breaker  coefficient,  and  wave-height  to 
water-depth  ratio.  Good  agreement  was  also  found  between  field  and 
laboratory  turbulence  intensities  when  normalized  by  the  bore  model.  The 
large  ratios  of  field-to-laboratory  wave  heights  (e.g.,  30:1)  and  frequencies  (e.g., 
10:1)  indicate  that  bore-model  scaling  of  surf  zone  turbulence  intensity  is 
robust  over  large  differences  in  scale. 


6.  APPENDICES 


A.  Bubbles 

The  signal  from  a  constant-temperature  hotfilm  anemometer  is 
corrupted  by  bubbles  because  air  has  much  lower  heat  capacity  than  water. 

The  hotfilm  power  (and  therefore  output  voltage)  required  to  maintain 
constant  temperature  is  much  smaller  in  air  than  in  water  at  the  same 
velocity,  [Delhaye,  1969].  The  hotfilm  output  voltage  signal  during  the  passage 
of  an  air  bubble  is  characterized  by  a  steep-sided  trough  (Figure  A).  Resch, 
Leutheusser,  and  Alemu  [1974],  studying  turbulence  levels  in  a  hydraulic 
jump,  used  dropout  voltage,  the  distance  between  the  local  signal  maximum  in 
the  water  and  the  minimum  in  the  bubble,  to  recognize  bubbles.  However, 
there  is  no  consistent  dropout  voltage  in  the  surf  zone  due  to  widely  varying 
bubble  sizes  and  advective  velocities,  so  we  used  a  slope  threshold  method 
from  Wang  [1985].  The  maximum  slope  of  the  signal  dropout  is  steeper  than 
voltage  changes  associated  with  flow  fluctuations.  The  hotfilms  were  sampled 
at  2048  Hz  for  bubble  removal.  Once  recognized,  each  bubble  signal  was 
replaced  using  linear  interpolation  and  the  time  series  reduced  by  block¬ 
averaging  to  64  Hz. 

B.  Calibration 

The  electromagnetic  current  meters  (calibrated  in  a  laboratory  flow 
channel  with  accuracy  of  a  few  cm/s;  see  Guza,  Clifton,  and  Rezvani,  [1988]) 
were  used  to  calibrate  the  hotfilms  in  situ.  Both  the  hotfilm  voltage,  which  is 
inherently  rectified  because  the  probe  senses  only  the  velocity  magnitude,  and 
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the  absolute  value  of  the  EM  current  meter  cross-shore  velocity  were  reduced 
to  8  Hz.  To  reduce  temporal  lag  problems  caused  by  the  small  distance 
separating  the  hotfilm  and  EM  current  meter,  the  data  points  in  each  of  the 
two  time  series  was  independently  ranked  according  to  magnitude  and  then 
recoupled  so  that  each  hotfilm  voltage  was  paired  with  the  EM  speed  of  the 
same  rank.  The  data  pairs  were  then  averaged  for  bins  spanning  5  cm/s.  A 
log-log  polynomial  provided  an  acceptable  and  convenient  fit  to  the 
resulting  paired  data.  Each  512  s  hotfilm  time  series  was  independently 
calibrated  and  plotted  (Figure  B(l)).  Time  series  with  obviously  bad  fits  were 
discarded.  The 

EM  and  calibrated  hotfilm  time  series  were  similar  at  sea-swell  frequencies 
(Figure  B(2)). 

C.  Binomial  Expansion  of  Hotfilm  Signal 

A  hotfilm  probe  senses  and  rectifies  flow  from  all  directions,  giving  a 
signal  of 

where  u  =  U  +  u'.  If  U»u'  and  U»V,W  then 

5(0-((C/+ay+v^+w^)*/^ 

so  that 

Sit)  =  (f/2  +  2Uu'  +  Oiu'^))^^  ^  U  ^  u’  ^  Oiu'^). 
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7.  FIGURES 


Figure  1.  Instrument  package  mounting  frame  with  positions  of 
hotfilms  (HF),  electromagnetic  current  meters  (EM),  and  pressure 
sensor  (P).  Dimensions  are  in  cm.  (a)  Plan  view,  (b)  View 
looking  cross-shore,  (c)  Viev'  looking  alongshore. 
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Figure  3.  Representative  time  series:  (a-d)  In  surf  zone;  60%  of 
waves  broken;  h=178cm;  Y=H/h=0.33.  (a)  Sea  surface  elevation 
from  the  pressure  data  and  linear  theory,  (b)  Hotfilm  speed, 
(c&d)  Hotfilm  high-pass  filtered  at  f„=2.35  Hz,  (c)  Elevation  aoove 
bed:  1^=119  cm.  (d)  ^=84  cm.  (e-g)  Seaward  of  surf  zone;  3%  of 
waves  broken;  h=149  cm;  'y=:H/h=0.28;  {;=119  cm.  (e)  Sea  surface 
elevation,  (f)  Hotfilm  speed,  (g)  Hotfilm  high-pass  filtered  at  f;.. 
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Frequency 

Figure  4.  Frequency  spectrum  of  a  typical  512  s  hotfilm  time 
series,  showing  orbital  wave  (F'^  slope)  and  inertial  subrange 
turbulent  slope)  regimes.  h=134  cm;  ^=0.29;  DOF 
increases  with  frequency  from  62-11000. 
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Figure  5.  Dissipation  rates  (based  on  1/8  s  subsections  of  a  512 
s  record)  plotted  on  semilog  paper  vs  probability  (e^c).  The  right 
ordinate  is  stretched  so  that  data  from  a  lognormal  distribution 
lie  on  a  straight  line.  The  corresponding  "standardized  normal 
random  variable"  (e-p)/a,„£  is  shown  on  the  left  ordinate.  A 
straight  line  (ideal  lognormal  distribution)  is  fit  through  the  data. 
The  inverse  slope  of  the  line  is  and  the  zero-crossing  is  used 
to  estimate  p. 
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Figure  8.  Vertical  profiles  of  turbulence  intensity,  u',  scaled  by  (gh)  Field  mean  and 
bounds  (solid  linos)  are  from  f’igure  7(c).  Laboratory  results  are  from  (a)  Stive  1 19801  and 
(b)  Hattori  and  Aono  [1985].  The  laboratory  data  are  sorted  by  (a)  h/boi,,  the  ratio  of  local 
depth  to  breaking  depth  and  (b)  x/x,j,  the  ratio  of  distance  from  shore  to  surf  zone  width. 


s  from  (a)  Stive  1 1980]  data  aiul 
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Figure  B.  (1)  Typical  hotfilm  calibration  curve.  Error  bars  are  ±1 
standard  deviation  in  each  5  cm/s  bin.  (2)  Cross-shore  velocity 
from  the  EM  current  meter  and  hotfilm  (derectified  to  match  the 
sign  of  the  EM). 
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